The dynamics of a classical spinning particle in Vaidya space-time 
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Based on the Mathisson-Papapetrou-Dixon (MPD) equations and the Vaidya metric, the motion 
of a spinning point particle orbiting a non-rotating star while undergoing radiation-induced gravi- 
tational collapse is studied in detail. A comprehensive analysis of the orbital dynamics is performed 
assuming distinct central mass functions which satisfy the weak energy condition, in order to de- 
termine a correspondence between the choice of mass function and the spinning particle's orbital 
response, as reflected in the gravitational waves emitted by the particle. The analysis presented here 
is likely most beneflcial for the observation of rotating solar mass black holes or neutron stars in 
orbit around intermediate-sized Schwarzschild black holes undergoing radiation collapse. The possi- 
bility of detecting the effects of realistic mass accretion based on this approach is considered. While 
it seems unlikely to observe such effects based on present technology, they may perhaps become 
observable with the advent of future detectors. 

PACS numbers: 04.30.-w, 04.30.Db, 04.25.-g, 04.70. Bw 



I. INTRODUCTION 



One of the earliest predictions of Einstein's theory of 
gravitation is the existence of gravitational waves due to 
the motion of massive sources. While strong indirect ev- 
idence [1] exists via observations of the Hulse- Taylor bi- 
nary pulsar PSR 1913-|-16-whose orbital period decreases 
in agreement with general relativity's predictions when 
accounting for the emission of gravitational radiation- 
direct observation of the phenomenon has not yet been 
achieved. If gravitational waves exist in the manner pre- 
dicted by general relativity, then by a detailed analysis of 
the shape and amplitude of the generated waveforms, it is 
theoretically possible to extract useful information about 
the physical properties of the collective mass source that 
generated the radiation. 

The search for gravitational waves is a topic of in- 
tense activity, and already a number of detectors have 
been built or are under development for the purpose of 
observing gravitational radiation. Arguably the most 
prominent of the gravitational wave detectors in exis- 
tence or under development are the ground-based Laser 
Interferometer Gravitational Wave Observatory (LIGO) 
[2, 3] and the space-based Laser Interferometric Space 
Antenna (LISA) [4]. These two facilities complement 
each other in that LIGO is sensitive to high-frequency 
gravitational waves of 10 — 1000 Hz generated from the 
inspiral of equal-mass binary neutron stars or black holes, 
while LISA is designed to detect lower frequency waves 
of 10~* — 10~^ Hz corresponding to the motion of solar 
mass neutron stars or black holes in orbit around super- 
massive black holes. Other ground-based interferometers 
that exist worldwide are GEO 600 [5], TAMA 300 [6], and 
VIRGO [7] . While this form of radiation is notorious for 
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being extremely difficult to detect, gravitational waves 
offer a valuable opportunity to provide insight about the 
sources which generate them. This includes the earliest 
moments in the known Universe that are not accessible 
by other observational means, since they can propagate 
over literally cosmological distances with no noticeable 
dispersion whatsoever [4]. Gravitational radiation then 
offers the possibility of gleaning observational informa- 
tion about the components that drive early Universe cos- 
mology within the first second of the Big Bang, such as 
evidence for string interactions and cosmic inflation. In 
contrast, information gained from electromagnetic radi- 
ation is confined to approximately 300,000 years after 
the Big Bang, when temperature fluctuations in the cos- 
mic microwave background radiation can only provide, 
at best, some insight about the formation of large-scale 
structures in the known Universe, with no solid informa- 
tion about the dynamical processes before the surface of 
last scattering. 

An especially interesting area of general relativity re- 
search concerns how ordinary stars of sufficient mass 
undergo gravitational collapse to form neutron stars or 
black holes. Starting from the first quantitative "dust" 
model introduced by Oppenheimer and Snyder [8] , which 
assumes a perfect fluid with zero pressure, more sophis- 
ticated models [9] now exist, each with their own unique 
features that depend on the star's internal properties. 
It appears that gravitational collapse can be categorized 
into two main types. The first type concerns stellar 
core collapse where the central mass remains constant, 
while the second type involves the accretion of mass onto 
the star. This second form of gravitational collapse has 
the particularly interesting effect of making the exterior 
space-time background dynamic with the increase of stel- 
lar mass, with possible consequences for the motion of 
any object in its orbit. Under suitable conditions, such 
an object may sense a noticeable distortion of the space- 
time background that can result in some nontrivial re- 
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sponse from an otherwise stable elliptical orbit. If the 
orbiting object has sufficient mass, then the gravitational 
radiation it generates from this dynamical response may 
possibly become observable by present or near-future de- 
tectors under development. 

Another interesting area of general relativity research 
concerns the dynamics of macroscopic objects with spin 
angular momentum propagating in a curved space-time 
background. Because most astrophysical objects have 
acquired at least some spin angular momentum during 
their formation, it is very important to know how such 
objects interact in strong gravitational fields. This prob- 
lem was first examined many years ago by Mathisson [10] 
who attained a leading order correction to the geodesic 
equations in the form of a direct spin coupling to the 
Riemann curvature tensor. His work was reformulated 
by Papapetrou [11], who described the spinning particle 
in terms of a matter field confined within a thin world 
tube in space-time. Many approaches have been taken 
to describe the motion of extended objects with spin 
[12]. While they identify the leading order correction 
for the equations of motion, they differ in the contribu- 
tions due to higher-order multipole moment expansions 
of the spinning mass. Arguably the most conceptually 
successful model was put forward by Dixon [13], who pro- 
posed a self-consistent set of equations which incorporate 
to infinite order all multipole moment expansions of the 
spinning object. Regardless, for most practical purposes, 
truncating the equations of motion to leading order in the 
spin interaction is sufficient for a good basic description 
of the object's physical behaviour in strong gravitational 
fields, so long as the spinning particle is small in dimen- 
sion compared to the space-time's local radius of curva- 
ture. This is often referred to as the "pole-dipole approx- 
imation," first introduced by Mathisson and Papapetrou. 
A wide variety of research on classical spinning particles 
over the past few years have been conducted, which in- 
clude both formalism considerations [14], particle kine- 
matics and dynamics in generally curved space-time [15], 
and scattering in gravitational wave space-times [16]. 

Recent applications of these macroscopic spin models 
have led to some interesting consequences when applied 
to black hole space-times. One example involved study- 
ing the motion of the spinning particle's near a Kerr black 
hole along surfaces of constant azimuthal angle as defined 
by the space-time's symmetry axis [17]. A related exam- 
ple studied the gravitational waveforms generated by a 
spinning particle while falling into a Kerr black hole from 
infinity [18] and in circular orbit [19]. A particularly in- 
teresting example concerned the search for evidence of 
deterministic chaos within the gravitational waves gen- 
erated by a spinning point particle in orbit around a 
Schwarzschild black hole [20, 21]. Such a claim would 
have obvious implications for the ability of existing and 
future gravitational wave detectors to extract any mean- 
ingful information from the predicted waveforms. How- 
ever, a more recent study of the evidence for chaos in 
a Kerr background suggested that the spins required to 
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FIG. 1: Initial spin orientation {6, (p) for the spinning particle. 
For the purposes of this paper, the observer is located along 
the X-axis. 



obtain chaotic motion are not likely to be found for as- 
trophysically realistic systems [22, 23]. 

One possible application that has apparently not re- 
ceived any attention so far concerns spinning particle in- 
teraction with a star while undergoing gravitational col- 
lapse into a black hole, with the possibility of observing 
the phenomenon through gravitational radiation. This 
may perhaps be due to the relative unlikelihood of find- 
ing a candidate system that captures black hole forma- 
tion in progress, compared to fully formed black holes. 
Nonetheless, the usefulness of looking into this possibility 
seems clear, given the great precision required to gener- 
ate templates in order to compare with a possible signal. 
Because the spinning particle can act like a probe of the 
space-time background, the gravitational radiation from 
such a particle has the ability, in principle, to refiect the 
dynamical properties of the collapsing star while in orbit, 
and that the possibility exists to find signatures of this 
effect in the gravitational radiation emitted from it. 

The purpose of this paper is to examine the dynami- 
cal motion of a spinning point particle described by the 
Mathisson-Papapetrou-Dixon (MPD) equations in orbit 
around a non-rotating compact object under gravita- 
tional collapse due to radiation accretion. The space-time 
background is then described by the Vaidya metric [24] 
in terms of an advanced null co-ordinate to account for 
the infall of radiation onto the star's surface. Although 
the metric is very simple in form, it nonetheless has valu- 
able properties for cleanly identifying where the precise 
effects of gravitational collapse can appear in the -I- and 
X polarization modes, particularly when the point par- 
ticle's spin interaction with the background is taken into 
account. This paper begins in Section II with the MPD 
equations of motion for a spinning point particle, followed 
in Section III by a detailed numerical analysis of the 
particle's orbit and gravitational radiation output in the 
Vaidya background, assuming various mass functions de- 
fined by monotonically increasing elementary functions. 
The underlying purpose of Section IV is to establish the 
response pattern between the mass function's growth rate 
and the orbiting particle's motion, particularly when re- 
flected in the gravitational waveforms. This is followed 
by Section V, which examines the observational possibil- 
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FIG. 2: Central mass functions considered for this paper, 
where e = 0.01. 



ities that follow from this analysis, particularly concern- 
ing rotating solar mass black holes and neutron stars in 
orbit around intermediate-sized black holes. The likeli- 
hood of detecting the effects of mass accretion based on 
realistic astrophysical models is also considered here. A 
brief conclusion in Section VI then completes this paper. 

The space-time metric is described in +2 signature, 
where calculations follow the conventions adopted in the 
reference by Misner, Thorne, and Wheeler [25], and geo- 
metric units of G = c = 1 are assumed throughout. 



II. EQUATIONS OF MOTION FOR A 
SPINNING POINT PARTICLE 



For a spinning test particle, the MPD equations of mo- 
tion are 



dr 
dr 



1 



V" S' 



(la) 
(lb) 



where = dx^ /dr is the particle's centre-of-mass four- 
velocity defined with respect to some affine parameter r 
that is not necessarily proper time, P^ is the particle's 
four-momentum, and S''"' is the particle's antisymmetric 
spin tensor. The tensor components are defined with 
respect to a co-ordinate basis. On their own, (la)-(lb) do 
not completely determine the particle's motion, as they 
require two constraints to be satisfied. The first type 
is an equation to specify the orthogonality conditions of 
the spin tensor. Following Dixon, the choice adopted 
for the spin condition is S*"^ Pp = Q, which allows for 
explicit expression of V^^ in terms of and S^'^ [26]. 
This results in 
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where m ~ (— P^P^) ' is identified as the point 
particle's mass. The second type of constraint is a 
parametrization condition required to determine P • V 
in (2), where the choice adopted here is to set V'^ = 



dt/dr = 1 throughout the particle's motion. This al- 
lows for identification of the trajectory's time develop- 
ment with the space-time background's co-ordinate time. 

The space-time background used to describe gravita- 
tional collapse from infalling null shells is the Vaidya met- 
ric [24] 



ds^ = - 1- 



2M{v) 



dv^ + 2dvdr + dn^ (3) 



written in Eddington-Finkelstcin form, where dOF' = 
dO^ + sin^ 9 dip^ and u — j/(t, r) is the advanced null 
co-ordinate. Given co-ordinates X^' = [v, r, 6*, ip) cor- 
responding to (3), it follows from the Einstein field equa- 
tions that Rnu — Stt Tan for a null fluid source, where 
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dv 
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and = (1, 0, 0, 0) is the ingoing null generator in these 
co-ordinates satisfying P^ = 0. 

While it is possible to directly use (3) for the particle 
motion's time development with respect to v, the metric 
can be expressed in terms of ~ {t, r, 6, ip) co-ordinates 
in order to make use of the co-ordinate time t. To do this 
involves using the tortoise co-ordinate condition 
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to write dv as a linear combination of dt and dr, where 
Mq is the static Schwarzschild mass of the central body. 



Then by the co-ordinate transformation 
(6), the Vaidya metric takes the form 
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dt dr + r^ dn^ (7) 



in terms of lapse and shift vector functions, with the 
attractive feature that (7) reduces to the Schwarzschild 
metric for M{v) Mq. The new ingoing null generator 
in (i, r, 6, ip) co-ordinates becomes 
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and the corresponding Ricci tensor for (7) is then 

R^u — l^^fi^i^: (9) 

0. 



where it can be confirmed that 1^- 
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FIG. 3: Orbital trajectories for a spin- 
ning particle around a Schwarzschild 
black hole. Fig. 3(a) describes the 
orbit for a particle without MPD 
spin-gravity coupling (A = 0), while 
Fig. 3(b) has the coupling included 
(A = 1). The initial spin orientation 
is 6 = if — tt/A for both cases. 



(a) A = 



(b) A = 1 






(a) log (A = 0) 



(b) tanh (A = 0) 



(c) exp (A = 0) 






(d) log (A = 1) 



(e) tanh (A = 1) 



(f) exp (A = 1) 



FIG. 4: Orbital trajectories for a spinning particle around accreting masses described by different mass functions. For all cases 
considered here, e = 0.10 and the initial spin orientation \& Q = (p = tt/A. Figs. 4(d)-4(f) indicate that the particle's polar 
angular position oscillates about 9 — n/2 due to MPD spin-gravity coupling, as observed in Fig. 3(b). 



III. FRAMEWORK FOR NUMERICAL 
SIMULATION 

A. Equations of motion 

Before proceeding to any numerical analysis of the 
MPD equations (la)-(lb), it is useful to express them 
in a purely vector form, since issues concerning round- 
off errors become more manageable than otherwise. This 
requires writing the equations in terms of a spin vector 
S", defined as 



1 
2^ 



(10) 



where = \/~9 (^tivap is the Levi-Civita tensor de- 

scribed in terms of the alternating symbol cr^^ap with 
o'oi23 = 1- Then the spin tensor in terms of 5" is 

S"'^ = —t°'P^"'P^S^, (11) 
m 

where e'^''"^ = — -y= a^^"'^ . The spin condition 

S'^P^ = {) (12) 

then naturally follows from (10), while the spin magni- 
tude S is defined by 



S"'^ Sap, 



(13) 
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with units of mA4{v). 

Besides the MPD equations, the time evolution of v 
must also be taken into account, since information about 
the mass function M{v) has to be related to the par- 
ticle's position with respect to co-ordinate time t in a 
consistent way. This requires an extension of the full 
set of equations to incorporate the algebraic constraint 
given by the tortoise co-ordinate condition (6), resulting 
in a Type I differential-algebraic system of equations to 
solve [27]. Therefore, the extended MPD equations for 
numerical integration are 
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(14c) 
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where the spin condition (12) is assumed throughout the 
particle's motion, (6) is utilized to describe (14d), and A 
is a dimensionless parameter added to tag all terms with 
MPD spin-gravity coupling. This turns out to be a very 
useful device, since the tuning of A allows for direct com- 
parison between spin evolution involving strictly parallel 
transport of the spin vector (A = 0) with the contribution 
due to spin-gravity coupling (A = 1). 



mass functions arc predefined to remain constant from 
the start of the simulation until v — vq, when the mass 
accretion rate is then "switched on" and M{v) is allowed 
to increase by a small relative value e < 1 over some 
predetermined lifetime tq before it gets "switched off". 
While there is no claim that this construction for the 
mass functions accurately reflects physically realistic sce- 
narios, the point of this exercise is to demonstrate the 
robustness of the MPD equations to describe the effects 
of the particle's spin on its orbital dynamics, depending 
on the growth rate of the mass increase. 

All dimensional parameters are described in terms of 
Mq, the mass of the central body prior to the start of 
mass accretion at v = vq. In units of Mq, the mass 
functions considered are the following: 

M{v) = 1, (15a) 
= l + £ln[l + (j.-j.o)/ro], (15b) 
= l + e [exp [{v - vq)/tq\ - 1] , (15c) 

= 1 + 1 {tanh [{u - vo)/to] + 1} , (15d) 

where (15a) corresponds to a Schwarzschild black hole, 
(15b) describes a logarithmic mass increase, (15c) refers 
to an exponential increase, and (15d) defines an increase 
due to a hyperbolic tangent function. For future refer- 
ence, these functions are respectively denoted by sch, 
log, exp, and tanh in the plots presented in this pa- 
per. All simulations presented for these functions assume 
To = 10^ Mq, and the choice of vq for each mass function 
is such that M{u) = 1 + e/2 at around t = 10"* Mq. The 
relatively large choice for tq is to allow for the mass func- 
tions to gradually increase over several orbital cycles. In 
addition to these mass functions, a step function mass 
increase (denoted hereafter by step) is simulated by us- 
ing (15d) with To = IMq. a plot of the mass functions 
over time is presented in Figure 2. 



B. Numerical details 

The integration is performed in Fortran 90, using 
a Runge-Kutta fourth-order algorithm with fixed time 
steps. For each simulation nm presented here, the par- 
ticle is initially positioned in the x-direction of the co- 
ordinate frame defined by the central mass function lo- 
cated at the origin, and follows a counterclockwise orbit. 
This is illustrated in Figure 1, where the particle's initial 
spin orientation is defined by {0. 0) co-ordinates follow- 
ing the standard definition of spherical co-ordinates with 
respect to the particle's local frame. 

The central mass for the Vaidya metric can be de- 
scribed by an arbitrary fmiction of v, so long as it is 
monotonically increasing to preserve the weak energy 
condition. With this in mind, a variety of mass func- 
tions are considered in the hope of establishing a corre- 
spondence between the time development of M{u) and 
the dynamical response of the spinning point particle. 
For the purpose of making well-defined comparisons, the 



IV. NUMERICAL ANALYSIS 

With the equations of motion (14a)-(14d) and mass 
functions (15a)- (15d), it is possible to determine a wide 
variety of kinematic and dynamical properties associated 
with the spinning point particle. This section details the 
numerical results that are obtained from this analysis. 
All simulations considered here assume for initial condi- 
tions a particle mass of m = IO^^Mq. This is consis- 
tent with the assumption that the spinning particle is 
very small compared to the central body in order to sat- 
isfy the test particle assumption. The spin-to-mass ratio 
considered here, which also satisfies the required approx- 
imation S/(mr) <C 1, is S/m = 5 x 10~^ Mq, with ex- 
ceptions noted as required. As with the choice of m, 
the spin magnitude is sufficiently large to accommodate 
for realistic astrophysical conditions while avoiding nu- 
merical instabilities that follow from making S too large. 
Except where otherwise indicated, e = 0.01 and the ini- 
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(b) log 



(c) step 



FIG. 5: Particle displacement oflt the orbital plane for Schwarzschild, logarithmic, and step mass functions, where e = 0.01. 
The dotted line corresponding to the largest amplitude refers to 9 = n /2 and ip — 0, the next largest amplitude is described by 
the thin solid line for 6 — ip = 7r/2, and the thick solid horizontal line corresponds to = 0. 
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FIG. 6: Effects of initial spin orientation on the particle's or- 
bital stability, assuming the hyperbolic tangent mass function 
with £ = 0.01, V = 0.395 c, and 5/m = 5 x 10"^ Mo. The solid 
line corresponds to A = 0, while the dotted line refers to A = 1. 
From Fig. 6(a), it is clear that spin orientation along ^ = 
leads to orbital stability when MPD spin-gravity coupling is 
taken into account, while from Fig. 6(b) an orientation along 
9 — TT leads to an accelerated orbital inspiral. 



tial orbital speed of the particle is i; = 0.350 c, which 
corresponds to an initial radius of r « 8 Mq- 



A. Orbital kinematics 

A sampling of orbital trajectories for a classical spin- 
ning particle are presented in Figures 3 and 4, where 
Figure 3 describes the orbit around a Schwarzschild 
black hole and Figure 4 lists the orbits for the logarith- 
mic, hyperbolic tangent, and exponential mass functions, 
where e — 0.10. As well, the initial spin orientation is 
9 — ip ~ tt/A for each plot. 

From Figure 3(a), it is clear that the particle follows 
a stable circular orbit due to purely geodesic motion 
(A = 0), while Figure 3(b) shows the contribution of 
the MPD spin-gravity coupling (A = 1) with a slight 
quasiperiodic oscillation of the particle's position about 
the 9 = 7t/2 orbital plane. However, for these initial 
conditions the orbital radius remains largely unchanged 
due to the force and torque terms in the MPD equations. 
In contrast, Figure 4 shows that the accreting mass can 
significantly alter the trajectory for sufficiently large e, 
leading to precessing ellipsoidal orbits. In particular, it is 
clear from comparing Figures 4(a)-4(c) to Figures 4(d)- 
4(f) that the influence of MPD spin-gravity coupling is 
relatively small compared to the particle's centre-of-mass 
motion due to the geodesic equation. Nonetheless, there 
are differences in that there is both some oscillation about 
the orbital plane and a slight reduction of the orbital ec- 
centricity due to the presence of spin. Not surprisingly, 
the orbiting particle responds in proportion to the degree 
of central mass increase with time, and that the MPD 
force and torque terms arc similarly responsive for the 
assumed spin orientation. 

To determine the effect of MPD spin-gravity coupling 
on particle confinement to the orbital plane, Figure 5 
describes its z-component amplitude over time for the 
Schwarzschild, logarithmic, and step mass functions. For 
spin orientation prepared along the radial direction {9 = 
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FIG. 7: Polar angular displacement for Schwarzschild, log- 
arithmic, and step mass functions, where e — 0.01 and 
9 = (f = n /A. The oscillating solid and dashed lines cor- 
respond to A = 1, where the latter describes the displacement 
due to a Schwarzschild black hole. In contrast, the solid hor- 
izontal line for A = shows that MPD spin-gravity coupling 
is the sole contributor to the polar displacement observed. 
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FIG. 8: Azimuthal angular displacement for Schwarzschild 
and hyperbolic tangent mass functions, where e = 0.01. For 
Fig. 8(b), the effect of MPD spin-gravity coupling leads to a 
decrease in angular displacement compared to geodesic mo- 
tion (A — 0) due to 9 = spin orientation, as described by the 
dotted line. The dashed line corresponding to ^ = tt shows a 
slight increase by roughly the same amount. 



7r/2 and ip = 0), the amplitude is z 0.3 Afg, while an 
orientation along the tangential direction [6 = (p — 7r/2) 
produces a slightly reduced amplitude of z « 0.2 Mq. 
This leads to a 2.50-3.75% displacement when compared 
to the particle's initial radial position. It is also clear 
from Figures 5(a)-5(c) that a spin orientation normal to 
the orbital plane {6 ~ 0) produces no obvious displace- 
ment. In addition, there is no evidence of amplitude 
growth due to mass accretion, regardless of the particle's 
spin orientation. 

It is of interest to examine the influence of spin on the 
particle's orbital stability over time for situations where 
the centre-of-mass motion alone leads to an inspiral. This 
is demonstrated in Figure 6 for the hyperbolic tangent 
mass function, which compares the lifetime of the parti- 
cle's orbit, where v = 0.395 c. Comparison of Figures 6(a) 
and 6(b) shows that the spin component along the 9 = 
direction is solely responsible for producing stability in 
the particle's orbit, and that a spin component directed 
along the 9 = n direction increases the inspiral rate when 



compared to purely geodesic motion. Given that most (if 
not all) observed binary systems have their orbital and 
spin angular momenta pointed in the same general direc- 
tion, this result makes intuitive sense. 

As demonstrated by the orbital trajectories of Fig- 
ures 3 and 4, there is an oscillation of the particle's 
polar angular position about 9 — t:/2 due to the force 
and torque terms in the MPD equations of motion. A 
closer examination of this effect is presented in Figure 7, 
with a comparison of the polar angle displacement for 
Schwarzschild, logarithmic, and step mass functions as 
a function of time, with e — 0.01. All plots here show 
that a polar displacement of A9 w 7.7 x lO"'^ tt about 
the 9 = 7r/2 orbital plane occurs due to MPD spin- 
gravity coupling. When mass accretion is taken into 
account, the frequency of oscillation increases in accor- 
dance with the rate of mass increase, as suggested by 
Figures 7(a) and 7(b), though the maximum displace- 
ment from 9 = 7r/2 remains largely constant. 

It is also evident from Figure 8 that the MPD spin- 
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FIG. 9: Phase difference in particle's centre-of-mass azimuthal position due to dynamical central mass when compared to 
Schwarzschild mass for A = 1. Figs. 9(a)-9(c) show that mass accretion induces a phase advance in the particle's orbit that 
reflects the rate of mass increase. 
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FIG. 10: Polar component of the particle's four-momentum 
for the logarithmic and hyperbolic tangent mass functions, 
where e — 0.01. The dashed line describes the amplitude due 
to a Schwarzschild black hole. 



gravity coupling has an influence on the overaU angular 
displacement in the azimuthal direction, depending on 
the initial spin orientation. For spin orientation with a 
component along 9 = 0, the particle's orbital frequency 
decreases when compared to that due to centre-of-mass 
motion only, while a component along 9 = tt drives up 



the frequency by about the same rate. This is demon- 
strated by Figure 8(b) for the hyperbolic tangent mass 
function. There is no observed difference in the particle's 
azimuthal position due to MPD spin-gravity coupling for 
spin orientations confined to the orbital plane. 

Another interesting feature is presented in Figure 9, 
which plots the degree of phase advance defined by 

Sip{t) EE </?(t)vaidya " <i2(i)sch (16) 

for a given moment in time during the simulation. Com- 
paring Figures 9(a)-9(c) for the hyperbolic tangent, ex- 
ponential, and step mass functions, respectively, it ap- 
pears that the amount of phase advance over time is 
roughly proportional to the rate of mass increase at a 
given moment. Given that the orbital period is around 
146.5 Afo, a n/2 phase shift beyond what is plotted at 
< = 10^ AIo for the hyperbolic tangent function occurs 
after about 14 particle orbits, while the same effect oc- 
curs after roughly 10 orbits for the exponential function. 
Reflecting the sudden mass increase at t = 10'' Mq for 
the step mass function, the 7r/2 shift occurs in just flve 
orbits of the particle. 

B. Particle dynamics 

The extended MPD equations of motion (14) form a 
twelve-dimensional dynamical system for the time evo- 
lution of the particle's position, linear momentum, and 
spin angular momentum vectors. It follows that a com- 
plex interplay should exist between the linear and spin 
angular momentum components that leads to the orbital 
kinematics described earlier. While it is difficult to deter- 
mine precisely how the interplay unfolds, it is possible to 
see some interesting properties of the dynamical variables 
in isolation, particularly when there is a time-evolving 
metric due to mass accretion. 

One example is found in Figure 10, which describes 
the polar component of the particle's four-momentum 
over time. When mass accretion is taken into account, 
as shown by Figures 10(a)-10(b), the polar momentum 
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00> 



(a) sch (e = 0) 



(b) tanh (0 = 0) 



(c) tanh (e = tt) 



FIG. 11: Azimuthal component of the particle's four-momentum for Schwarzschild and hyperbolic tangent mass functions, 
where e — 0.01. Fig. 11(b) shows that the MPD spin-gravity coupling reduces the orbital speed for 6 — 0, while the opposite 
is true for ^ = tt, as shown by Fig. 11(c). 
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(a) sch (e = 0) 



(b) tanh {§ = 0) 



(c) tanh (e = tt) 



FIG. 12: Polar component of the particle's spin angular momentum for Schwarzschild and hyperbolic tangent mass functions, 
where e = 0.01. The mass accretion moderately boosts the magnitude of 5*", as shown by comparing Figs. 12(a)-12(b). The 
MPD spin-gravity coupling reduces the spin magnitude for 9 = 0, while Fig. 12(c) shows an increase in the magnitude for 



component oscillates more rapidly in response to where 
the mass accretion rate is strongest. A second example of 
note shows the time-averaged evolution of the particle's 
azimuthal component of its four-momentum, described 
by Figure 11. (All time-averaged plots presented in this 
paper are calculated over an interval of At = lO^Afo.) 
From Figures 11(b) and 11(c), it is evident that the MPD 
spin-gravity coupling acts to reduce the tangential mo- 
mentum if the spin component is projected along = 0, 
while it increases {P'^) if the spin has a projection along 
9 = IT. Similarly, from Figure 12, the polar component of 
the particle's spin angular momentum exhibits a similar 
behaviour. Comparing Figures 12(b) and 12(c), it follows 
that the MPD spin-gravity coupling reduces the magni- 
tude of (S^) for projections along ^ = while it increases 
the magnitude for projections of the spin orientation in 
the opposite direction. The effect of the mass accretion, 
as especially shown in Figures 11 and 12, is to moderately 
increase the overall magnitude of the particle's linear and 
spin angular momenta. 

It is particularly interesting to investigate the effects of 
mass accretion on the orbital angular momentum's rate 
of change with respect to time. A straightforward calcu- 



lation shows that the spinning particle's orbital angular 
momentum components are 

L'' = 0, (17a) 
= -r sin 61 P^, (17b) 
L'P = rP'^/ sine, (17c) 

while the components of the torque on the particle are 
"dt 



sin ( 





dt 



dt 



(18a) 



dl^ 
dt 



dL^ 
IT 



dP^ 

~dr 



2sme — P^ 
dt 



cos9{—P^ 
dt 



dip 



pi 



1 



sm( 



dP^ 
dt 



sm ti cos 



0^p. 
dt 



2^P' 
dt 



(18b) 



(18c) 



For a particle in circular motion confined to the 9 = tt/2 
orbital plane, (18) is identically zero, as expected. Fig- 
ure 13 is a time-averaged plot of (18a)-(18c) with respect 
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(a) sch (dV Idt) 



(b) sch {dL^ /dt) 



(c) sch (dL-^/dt) 





(d) exp (dV Idt) 



(e) exp (dL<^ /dt) 



(f) exp {dL'P/dt) 



FIG. 13: Time- averaged orbital torque on the spinning particle for Schwarzschild and exponential mass functions, where 
e = 0.01 and 6 = (p = it/4:. While some variation exists due to MPD spin-gravity coupling, Figs. 13(a)-13(f) show that orbital 
angular momentum is still conserved. 





FIG. 14: Radiation density for an 
observer comoving with the spinning 
particle due to the hyperbolic tangent 
mass function, where e = 0.01. 



(a) tanh {0 = 0) 



(b) tanh {§ = tt) 



to time for the Schwarzschild and exponential mass func- 
tions, where — ip — n/A. While there is some variation 
present in {dU /dt) due to MPD spin-gravity coupling, 
it is quite clear that orbital angular momentum is still a 
conserved quantity on average. 



C. Radiation density of the Vaidya Space-time 
background 

Because the Vaidya metric satisfies the Einstein field 
equations in matter described by a null fluid, a particle 
in orbit will come in contact with the surrounding radi- 
ation, with possible dynamical consequences. The radia- 
tion density corresponding to the Vaidya background, as 



measured by a local observer [24], is 



q = T^.V^V = ^R^^V^'V'', (19) 

O TT 



where R^i, is described by (9). 

It is of interest to find out whether the MPD equa- 
tions have any effect on the radiation density, as seen by 
a local stationary observer. To find out, a plot of (19) 
is presented in Figure 14 for the case of the hyperbolic 
tangent mass function. From comparing Figures 14(a) 
and 14(b) for initial spin orientations of ^ = and 9 = tt, 
respectively, it is shown that the MPD spin-gravity cou- 
pling induces about a 3% shift in the measured radiation 
density with respect to the maximum value determined 
by geodesic motion. It happens that only the compo- 
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nent of spin projected along the z-direction is sensitive 
to the difference in measured radiation density, and the 
results shown here are consistent with the kinematic ef- 
fects of the spin orientation presented earlier. Because of 
the velocity dependence inherent in (19), the local radi- 
ation density is analogous to a velocity-dependent drag 
effect with the dimensions of pressure. However, it is 
also apparent that with a maximum magnitude of only 
(J « 6 X 10""'^'' A^o~^ ^^^^ rapidly decays to zero, the radi- 
ation density of the Vaidya background appears to have 
little impact on the spinning particle's orbital motion. 



D. Gravitational waveforms and radiation 

The next step in this Section is to obtain the gravi- 
tational waveforms and generate the radiation output of 
the spinning particle while in orbit around the accret- 
ing central mass. For this exercise, the -|- and x polar- 
ization modes are described by the transverse traceless 
(TT) gauge, where the reduced quadrupole moment for- 
mula is assumed throughout for all relevant calculations. 
Although this assumption is most likely not sufficient for 
the purposes of comparing with potential observations 
by LIGO or LISA, where a more sophisticated analysis in 
terms of black hole perturbation theory may be required, 
it is still possible to gain valuable qualitative information 
about the waveforms' leading-order properties. 

To begin, the reduced quadrupole moment [21] for the 
point particle is 



(20) 



where Xj is the particle's Cartesian co-ordinate in space 
relative to the central mass. Before proceeding further, 
it should be noted that (20) is defined for co-ordinates 
that are isotropic in space, where the tensor indices are 
raised and lowered with the Minkowski metric. Because 
the Vaidya metric as represented by (7) is not isotropic 
in space, there will be some lack of precision in the grav- 
itational waveform amplitudes due to corrections in the 
radial and/or time co-ordinate to generate a more suit- 
ably defined metric. Using as a guideline the co-ordinate 
transformation [28] 



r(r) = 



r(r) 



r \ r J 



1/2 



(21a) 
(21b) 



to describe the Schwarzschild metric in terms of isotropic 
spherical co-ordinates {f,9,(p), it follows that f « 7 Mq 
for r « 8 Mo, with an error of approximately 12%. Be- 
cause it appears unlikely that (7) can be put into isotropic 
form in a simple manner, and given the likelihood that 
LIGO's or LISA'S signal-to-noise ratio will not notice the 



difference, the waveforms presented below offer a reason- 
ably suitable qualitative description of the gravitational 
radiation generated by the spinning particle. 

In the TT gauge and for an observer located at 
(ro, 6*0, </3o)) the + and x polarization mode amplitudes 
for the quadrupole moment formula are 



COS^ ^0 + 1 



[cos 2(po [hxx - hyy) 



1 



-I- 2 sin 2(^0 hxy] - - sin {Kx + hyy - h^z) 



- - sin 26*0 (cos (po hxz + sin <^o hyz) , (22a) 
hx = - i cos 6*0 [sin 2(pQ [hxx - Kv) ~ cos 2ipo hxy] 



+ sin 6*0 (sin (po hxz + cos (po hyz) , 



where 



hi. 



2 d Q ij 



(22b) 



(23) 



For the purpose of this paper, the observer is located 
along the orbital plane in the a;-direction, where 9o = 7r/2 
and (fo = 0. From this vantage point, any particle whose 
motion is confined to the orbital plane will only generate 
Therefore, any observation oi h^ with significant 
amplitude is a certain indicator of a force applied on the 
particle along the polar direction. 

To compare the contributions of MPD spin-gravity 
coupling to the -|- polarization mode, a set of plots for 
h+ for the Schwarzschild, exponential, and step mass 
functions arc listed in Figure 15, where 6 = if = n/A 
and £ = 0.01. A review of Figures 15(a) and 15(d) for 
the Schwarzschild mass shows that h+ reveals a slightly 
jagged peak structure in the waveforms when the force 
and torque terms in the MPD equations of motion are 
incorporated, though there is no obvious increase in the 
waveform's amplitude or frequency of oscillation. 

When the effects of mass accretion are included, the 
particle's orbital response does get reflected in the -|- 
mode, albeit with a very slight increase in the h+ ampli- 
tude over time. Comparison of Figures 15(b) and 15(e) 
with their respective Schwarzschild counterparts shows 
how h+ responds to mass accretion. This is evident in the 
step mass function examples of Figures 15(c) and 15(f), 
where it is obvious that a sudden change of mass had 
occurred at t = lO^Afo- It is particularly interesting to 
note for the step function examples that the effect of mass 
accretion alone is capable of producing a jagged peak 
structure in h+ due to purely geodesic motion (A = 0) 
that is virtually indistinguishable from the corresponding 
plot when MPD spin-gravity coupling is incorporated. A 
short spike does appear in the amplitude of Figure 15(f) 
when the mass accretion rate is almost infinite. Although 
it is possibly a numerical artifact only, the amplitude 
spike may be a legitimate signal due to MPD spin-gravity 
coupling, since the step mass function as defined from 
(15d) is formally smooth. 
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(a) sch (A = 0) 



(b) exp (A = 0) 



(c) step (A = 0) 




(d) sch (A = 1) 



(e) exp (A = 1) 



(f) step (A = 1) 



FIG. 15: Waveforms in the + polarization mode for Schwarzschild, logarithmic, and step mass functions, where e = 0.01. 
The initial spin orientation for all cases is 6 = ip = 7r/4. The effects of MPD spin-gravity coupling shows that a distinctive 
modulation appears in the waveform peaks. 
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(a) log (A = 


0) 



(d) log (A = 1) 



(b) exp (A = 0) 



1 lAfcl 

(e) exp (A = 1) 



1 

(c) tanh (A = 0) 



(f) tanh (A = 1) 



FIG. 16: Waveforms in the + polarization mode for logarithmic and exponential mass functions, where e = 0.01, 6 = <p = n/A, 
V = 0.395 c, and S/m = 5 x 10~''Afo. A well-defined amplitude modulation due to MPD spin-gravity coupling is present in 
Figs.l6(d)-16(f) during the final stages of particle inspiral. 
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(a) sch 



(b) log 



"i V 



(c) step 



FIG. 17: Waveforms in the x polarization mode for Schwarzschild, logarithmic, and step mass functions, where e = 0.01. The 
dotted line representing the largest amplitude corresponds to S = 7r/2 and ifi — 0, the thin solid line for the second largest 
amplitude refers to 8 — ip — 'k/2, and the thick solid horizontal line represents 8 — 0. 
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(a) log 



-f 5f ill' ' 



'''fi'''"'i'irir)'i'iirii.' 



'ir ll' ir ii' 



(b) exp 



FIG. 18: Waveforms in the x po- 
larization mode for logarithmic and 
exponential mass functions, where 
£ = 0.01, 6 = 0^ 7r/4, V = 0.395 c, 
and S/m = 5 x 10"^ Mq. 



A more interesting set of plots is found in Figure 16, 
which captures the waveform of a particle in the final 
stages of its inspiralling orbit. For both the logarith- 
mic and exponential mass functions considered, where 
V = 0.395 c, S/m = 5 x IQ-^Mq, and 6 = >f = 7r/4, the 
familiar-looking "chirp" signal appears during the par- 
ticle's final moments before its demise. Comparison of 
Figures 16(d)- 16(f) with their respective plots for A = 
suggests the presence of a modulation in the amplitude 
before final inspiral due to MPD spin-gravity coupling. 
However, the degree of modulation is very slight for these 
examples, and it is unlikely to be able to identify the pres- 
ence of particle spin in given the precision required 
for LIGO or LISA to identify a candidate signal. 

In contrast, of potentially much greater relevance for 
this paper is the effect of MPD spin-gravity coupling 
on the /ix waveform and the correlation of amplitude 
with its spin orientation. Because the particle oscillates 
about 9 = 7r/2, some x polarization radiation must be 
emitted, though it is not obvious about the amplitude 
strength when compared to as seen by an observer 
on the orbital plane. Figure 17 lists the x mode wave- 
forms for the Schwarzschild, logarithmic, and step mass 
functions, where the dotted line corresponds to initial 
spin orientation — 71/2 and = 0, the thin solid line 
describes 6 = (f> — 7r/2, and the thick solid horizontal 
line represents 6* = 0. For the case of an inspiralling 
particle. Figure 18 list the plots of ft-x for the logarith- 



mic and exponential functions, where v = 0.395 c, and 
S/m 5 X 10~^Mo. It is interesting to note that the 
maximum amplitude of /ix is about 5% of the corre- 
sponding amplitude for , as described by Figure 15. 
While the projection of spin orientation along the radial 
direction of motion (0 = 7r/2 and (p = 0) contributes the 
largest amplitude, the contribution along the tangential 
direction {9 = (p = 7r/2) is comparably large, as it gen- 
erates about 55% of the maximum amplitude. For spin 
orientation along the positive normal direction [9 = 0), 
there is no x mode radiation emitted. These results sug- 
gests the interpretation that, over time, the particle is 
trying to align its spin in the positive z-dircction as a 
point of stable equilibrium, reinforcing the same claim 
made earlier based on other data. 

A view of Figure 17(a) describing for the 

Schwarzschild mass shows a well-defined and symmetrical 
modulation of the amplitude due to the MPD spin grav- 
ity coupling. The effect of mass accretion, as described by 
Figures 17(b) and 17(c), causes a slight distortion in the 
amplitude, which becomes more noticeable with a faster 
rate of mass accretion. More significantly, however, the 
wavetrain gets compressed over the time scale where the 
mass accretion is most pronounced, with a more narrow 
modulation of h+. This is especially clear from the step 
function plots in Figure 17(c) where the particle is able 
to fit in roughly an extra six or more orbital cycles over 
the same time scale as compared to Figure 17(a). This 
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(a) sch (0 = = 7r/4) (b) log (0 = = 7r/4) (c) tanh {§ = ip = tt/4) 

FIG. 19: Gravitational wave amplitude for Schwarzschild and other mass functions, averaged over all directions. As shown by 
Figs. 19(b) and 19(c), the amplitude increases in direct response to the mass function's rate of growth. The MPD spin-gravity 
coupling slightly reduces the amplitude compared to the respective value for geodesic motion for each case. 
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(a) sch (A = 0) 



(b) tanh (A = 0) 



(c) step (A = 0) 



,1 f I, I. 



(d) sch (A = 1) 



(e) tanh (A = 1) 



(f) step (A = 1) 



FIG. 20: Energy emission rate for Schwarzschild, hyperbolic tangent, and step mass functions, where e — 0.01 and 8 — ip — 7r/4. 
The MPD spin-gravity coupling has a modulating effect on the emission rate, though for most mass functions considered the 
net energy loss is effectively zero. The only exception is for the step mass function, where there is a non-zero emission at the 
moment when the mass suddenly increases. From comparing Fig. 20(f) to Fig. 20(c), the MPD spin-gravity coupling slightly 
increases the energy emission rate for this choice of spin orientation, though the effect is still negligibly small. 



phase shift is also present for the logarithmic function in 
Figure 17(b), albeit a slightly more modest effect in this 



The next step is to examine the impact of mass ac- 
cretion on the predicted energy and orbital angular mo- 
mentum emission rates for radiation generated by the 
spinning point particle. One area where the impact be- 
comes immediately apparent occurs when plotting the 
overall amplitude A of the space-averaged gravitational 



wave [21], where from (23), 



2 /d^Qy d^gy^ 



(24) 



Plots of (24) scaled with rg are listed in Figure 19, which 
describe the gravitational wave amplitude for each of the 
mass functions. From Figure 19(a), the MPD spin grav- 
ity coupling has the effect of slightly reducing the overall 
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amplitude compared to that from purely geodesic motion, 
while generating a relatively large variation. When mass 
accretion is incorporated, it is clear from Figures 19(b) 
and 19(c) that the amplitude responds almost immedi- 
ately and in proportion to the mass accretion rate. It 
is interesting to note that, from Figure 19(c), the ampli- 
tude variation for the A = plot is roughly a quarter as 
for the A = 1 plot, and the reason for this anomalous be- 
haviour is not entirely clear, especially when compared to 
Figure 19(b). One possible explanation is that the spin- 
gravity interaction due to parallel transport gradually 
grows over time for a smoothly evolving mass function, 
while this opportunity is not available for the other mass 
functions considered. Therefore, it may be a numerical 
byproduct only. 

It is of interest to determine how much energy and or- 
bital angular momentum gets radiated away by the spin- 
ning particle in orbit around an accreting mass. Based 
on the quadrupole moment formula, the energy emission 
rate and the Cartesian component of the orbital angular 
momentum emission rate are [21] 

/dE\ _ 1 U^Q,, d^Q^i\ 

where these expressions are time-averaged over several 
cycles. The sign conventions adopted for (25) and (26) 
are that positive-valued energy and angular momentum 
are gained by the system. 

Figure 20 lists a set of energy emission rate plots for the 
Schwarzschild, hyperbolic tangent, and step mass func- 
tions, with spin orientation 9 = (f = n/A. Comparing 
Figures 20(a) and 20(d) for the Schwarzschild mass shows 
that the the MFD spin-gravity coupling modulates the 
emission rate, but remains essentially zero- valued. Simi- 
larly, Figures 20(b) and 20(e) corresponding to the hyper- 
bolic tangent mass function shows that while a gradual 
mass accretion adds some structure to the plot, the net 
effect on the emission rate is negligible. The exception 
to this rule occurs when the mass accretion rate is very 
rapid, such as found in Figures 20(c) and 20(f) for the 
step mass function, where a sudden burst of radiation is 
emitted from the particle at the moment where the mass 
increase rate is highest before settling back to zero emis- 
sion. The effect of MFD spin-gravity coupling for this 
instance is to reduce the emission rate to a practically 
negligible magnitude. 

The presence of particle spin adds an interesting di- 
mension to the decay of orbital angular momentum from 
gravitational radiation. This is evident from examining 
Figure 21, which details the time- averaged emission rate 
along all three spatial directions due to the Schwarzschild 
and exponential mass functions, assuming s = 0.01 and 
6 = (p = 7r/4 for the initial spin orientation. Not sur- 
prisingly, the loss of orbital angular momentum is over- 
whelmingly directed along the ^-direction, which is nor- 



mal to the = 7r/2 orbital plane, where the emission is 
two orders of magnitude higher than for the contributions 
along the radial and azimuthal directions. (It is impor- 
tant to remember that, for a particle in counterclockwise 
orbit along the 9 = it/2 plane, the orbital angular mo- 
mentum is directed along the negaMve ^?-direction.) 

Focussing for the moment on Figures 21(a)-21(c) for 
the Schwarzschild case, the first interesting point is that 
the system gains orbital angular momentum along the 
radial and azimuthal directions due to MPD spin-gravity 
coupling, while there is no corresponding emission for 
purely geodesic motion. The A = 1 plot for Figure 21(a) 
shows that \ {dL'' /dt) re.d\ ~ 3.5 x lO-^Mo. For Fig- 
ure 21(c), \r sm0{dLf/dt)rad\ « 2.4 x 10"* Mo, roughly 
an order of magnitude larger. However, the bulk of the 
radiation emission is along the ^^-direction, with magni- 
tude \r{dLydt)rs,d\ « 3.1 X 10-*^'Afo from Figure 21(b). 
Over the course of one orbit, the estimated radiation loss 
is 

\ dt 

w 5.5 X 10"^ Mo, (27) 

where At ~ 1.4 x 10^ Mq for a single orbit. By compar- 
ison, = -4.16 x 10-3 Mo and 5^ = -6 x 10"'"^ Mq 
for the spinning particle, so (27) radiates away a mag- 
nitude of angular momentum comparable to that of the 
particle's spin angular momentum. 

Concerning the contribution of MPD spin-gravity cou- 
pling to the angular momentum emission rate, this can 
be estimated by comparing the A = and A = 1 plots 
in Figure 21(e), where it follows that the spin interaction 
reduces the emission rate by about 5%. In addition, it 
follows that the exponential mass accretion rate increases 
{dL^ /dt)iad by roughly 10% over the time interval where 
the growth rate is in effect. 

E. Power spectrum analysis 

An important aspect of this analysis concerns deter- 
mining the frequency distribution of the power output 
associated with the gravitational radiation from the spin- 
ning particle. It is particularly interesting to know the 
impact of MPD spin gravity coupling on the power spec- 
trum, and its relative contribution compared to that 
of purely geodesic motion. Figures 22-24 describe the 
power spectra for a spinning test particle in stable quasi- 
circular orbit for each of the mass functions considered, 
assuming a spin orientation oi 9 = (p = w / A, where the 
spectra are defined in units of m^. 

For the Schwarzschild mass function with A = 0, Fig- 
ure 22(a) shows a single resonance peak at w 400 M^"^ 
with a maximum power of dE/du> « 8.8 x 10"'^ w^. This 
peak, corresponding to the frequency of the orbit, is com- 
mon to all power spectrum plots presented here. By 
comparison as shown in Figure 22(b), the contribution 
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(a) sch {dU' /dt)^^^ 



(b) sch {dL'^/dt),,^ 



(c) sch {dL'^/dt),^i 




(d) exp {dL^ /dt),,c,d 



(e) exp {dL'>/dt),^i 



(f) exp {dL^/dt),^i 



FIG. 21; Orbital angular momentum emission rate for Schwarzschild and exponential mass functions, where e = 0.01 and 
9 = if) = 7r/4. The dashed oscillating lines in Figs. 21(d) and 21(f) correspond to radiated orbital angular momentum of the 
particle orbiting a Schwarzschild black hole. 



(a) sch (A = 0) 




(b) sch (A = 1) 



FIG. 22: Power spectra for the Schwarzschild mass function, where 9 = (p = tt/4. 



of MPD spin gravity coupling has the effect of generat- 
ing two well-defined but smaller peaks, with a frequency 
separation of Au) « 100 Mq^. While there is a slight 
reduction of maximum power at the principal frequency, 
there is no fundamental difference present between the 
two plots. 

The effect of mass accretion imposes distinctive struc- 
ture on the power spectrum plots. As shown in Fig- 
ure 23, the the power spectra for the logarithmic and ex- 
ponential mass functions exhibit similar properties, both 
for A = and A = 1, particularly in that the principal 



peak is slightly broader than found in the Schwarzschild 
case. This is not surprising given that the mass func- 
tions evolve similarly over time. Comparison of Fig- 
ures 23(a) and 23(b) for the logarithmic mass function 
shows that the MPD spin-gravity coupling also generates 
the two smaller peaks with the same frequency separation 
as found in Figure 22, except that the lowest frequency 
one at a; « 300 Mq^ is slightly smaller in height, and that 
both smaller peaks are slightly distorted. These proper- 
ties also present for the exponential mass function case 
in Figures 23(c) and 23(d). 
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(a) log (A = 0) 



u.- [A/,7'1 

(b) log (A = 1) 





(c) exp (A = 0) 



u.' [A/,r'] 
(d) exp (A = 1) 



FIG. 23: Power spectra for the logarithmic and exponential mass functions, where 6 = ip — n /A. 



The hyperbolic tangent and step mass functions in Fig- 
ure 24 generate a power spectra with much greater com- 
plexity than found with the other functions. This is most 
evident when the effects of MPD spin-gravity coupling 
are considered, since at least four or five smaller peaks 
can be easily identified in Figures 24(b) and 24(d). The 
only exception in similarity comes from comparing Fig- 
ures 24(a) and 24(c) for A = 0, where the former reflects 
a much greater influence of MPD spin gravity coupling 
on the corresponding power spectrum than the latter. 



For the case of an infalling spinning particle with 
§ ^ 7r/4, V = 0.395 c, and S/m = 5 x IQ-^Mq, the 
peak structure for the power spectrum changes consid- 
erably. Figures 25 and 26 show that only one resonance 
largely appears and is shifted towards the viscinity of 
uj « 200 Mq^ for all plots. The only exception to this 
pattern is Figure 26(b), whose resonance peak remains 
at around uj « 400 M^j"^. For all cases, however, the 
power output at resonance drops by about an order of 
magnitude compared to the corresponding sets for stable 
orbits. 



V. OBSERVATIONAL POSSIBILITIES 

The essential motivation for examining the dynamical 
behaviour of spinning particles in Vaidya space-time is 
to determine the predicted gravitational waveforms that 
may, in principle, become observable by detectors such as 
LIGO or LISA. However, the likelihood of actually ob- 
serving the effects of radiation-induced mass accretion in 
this manner, both on conceptual and practical grounds, 
seems negligible when compared to the behaviour deter- 
mined by a Schwarzschild black hole. The main reason 
for this is because the amount of mass accretion required 
over a reasonable observation time is many orders of mag- 
nitude larger than allowed for by the Eddington luminos- 
ity limit for physically realistic mass accretion models. 
Further to this point is that because the central mass is 
non-rotating, the system lacks the presence of spin angu- 
lar momentum to drive up the mass accretion rate to a 
more robust level. 

To better illustrate the difficulty of observing spher- 
ically symmetric mass accretion in the waveforms, the 
Eddington limit imposes a mass accretion rate of [29] 



Af 3 X 10^22^-1 



Mo 

Mr.. 



(28) 
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FIG. 24: Power spectra for hyperbolic tangent and step mass functions, where e — 0.01 and 9 = (p = 7r/4. 



in geometric units, where 7 « 0.1 is the energy release 
efficiency of the outgoing photon flux to counterbalance 
the radiation mass accretion. To correct for the condition 
that only a small fraction of the infalling photons get 
radiated away, the mass accretion rate is then dM/dt = 
(1 — 7) M . Therefore, the time required to observe a 
central mass increase of AM = e A/q is 



At 



AAf 



dM/dt 
~ 5.8 X 10'' e years 



3.7 X 10^0 eM, 







(29) 



To then detect within the gravitational waveforms a 1% 
increase in the central mass requires an observation time 
on the order of 10^ years. Similarly, a one-year observa- 
tion time can yield, at best, an observable mass increase 
of only AM « 10~® Mq, which must be negligibly small 
for the design capabilities of LIGO or LISA. It is uncer- 
tain at present whether such a small effect may become 
reasonably detectable with the advent of future gravita- 
tional wave observatories. Nevertheless, the usefulness 
of this investigation comes from providing a theoretical 
testing ground for these same mass accretion models by 
comparing the predicted waveforms with actual signals. 

Given the power spectrum plots of Figures 22-26, it 
is important to know the size of central mass that is 



required for the resonance to fall within the frequency 
detection bands of LIGO and LISA. This can be deter- 
mined by letting Mq = aM^ for some unspecified pa- 
rameter a and evaluating the resonance peak at around 
« 400 Mq"^ in terms of M~^ = 2.03 x 10^ Hz. It then 
follows that 



.12 X lO'^a"^ 



Hz. 



(30) 



Assuming that the spinning point particle is a solar 
mass Kerr black hole or rapidly rotating neutron star, 
a — lO'^ and w ~ 8.12 x lO'* Hz, which is about 
two orders of magnitude beyond LIGO's upper limit in 
its frequency band, but may perhaps be open for ob- 
servation by a next-generation interferometer. To fall 
well within LIGO's bandwidth then requires that a = 
10^ — 10^, which implies that the orbiting point parti- 
cle must be an intermediate-sized Kerr black hole with 
TO = lO^^ — 10^ Mq. In order to satisfy LISA's bandwidth, 
it is necessary that a = 10^-10^, with m = 10^-10^ Mq. 
On the surface, it seems implausible to expect to find in 
abundance binary systems composed of two supermas- 
sive black holes. Nonetheless, it may not be completely 
outside the realm of possibility, since there appears to be 
no obvious reason why supermassive black hole binaries 
should not be found within the known Universe. 
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FIG. 25: Power spectra for logarithmic and exponential mass functions determined by a spinning particle in an inspiraling 
orbit, where e = 0.01, ^ = = 7r/4, i; = 0.395 c, and 5/m = 5 x 10"^ Mq. 



VI. CONCLUSION 

This paper investigated the dynamics of a classical 
spinning point particle in orbit around a gravitationally 
collapsing central mass due to radiation accretion, as de- 
scribed by the Vaidya metric. It is demonstrated that 
the Mathisson-Papapetrou-Dixon equations have a dis- 
cernable influence on the motion of the spinning particle. 
A comprehensive numerical analysis of the particle's or- 
bital kinematics and dynamics is performed for various 
mass functions which satisfy the weak energy condition. 
In addition, a leading-order description of the gravita- 
tional waveforms generated by the spinning particle is 
presented, along with a power spectrum analysis. For an 
observer positioned along the orbital plane, it is shown 
that the x polarization mode captures most clearly the 
properties of MPD spin-gravity coupling, with an ampli- 
tude of around 5% compared to the polarization mode 
for the choice of initial conditions considered. However, it 
is also shown that any effects for a realistic mass function 
with an accretion rate satisfying the Eddington luminos- 
ity limit must be very small. 

There are a number of potential considerations for fu- 
ture work based on this paper. Because the waveforms 



are determined using the leading-order quadrupole mo- 
ment formula only, it may be necessary to calculate more 
precisely defined waveforms based on black hole pertur- 
bation theory. It should be noted, though, that such an 
endeavour may yield only a minimum level of correction, 
given that a realistic mass accretion rate is very small 
in the first place. Likely a more meaningful advance of 
this work is to allow for rotation of the central mass, 
as described by the Kerr- Vaidya metric [30], and repeat 
the analysis presented here. Another interesting possible 
future development is to consider the spinning particle 
dynamics around a star undergoing a supernova explo- 
sion. While there is some complexity involved given that 
a supernova evolves in a highly asymmetric fashion, the 
fact that a sudden change of space-time curvature that 
follows from the event should have a significant impact on 
the orbital properties of the spinning particle, with con- 
sequences for the gravitational radiation emitted. These 
and other possible developments may be considered in 
the future. 
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FIG. 26: Power spectra for the hyperbolic tangent and step mass functions determined by a spinning particle in an inspiraling 
orbit, where £ = 0.01, § = tp = n/A, v = 0.395 c, and 5/m = 5 x 10"^ Mq. 
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